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ABSTRACT 

We analyze a suite of 33 cosmological simulations of the evolution of Milky Way-mass galaxies in 
low-density environments. Our sample spans a broad range of Hubble types at z = 0, from nearly 
bulgeless disks to bulge-dominated galaxies. Despite the fact that a large fraction of the bulge is 
typically in place by z ~ 1, we find no significant correlation between the morphology at z = 1 and at 
z = 0. The z = 1 progenitors of disk galaxies span a range of morphologies, including smooth disks, 
unstable disks, interacting galaxies and bulge-dominated systems. By z ~ 0.5, spiral arms and bars 
are largely in place and the progenitor morphology is correlated with the final morphology. We next 
focus on late-type galaxies with a bulge-to- total ratio B/T< 0.3 at z = 0. These show a correlation 
between B/T at z — and the mass ratio of the largest merger at z < 2, as well as with the gas 
accretion rate at z > 1. We find that the galaxies with the lowest B/T tend to have a quiet baryon 
input history, with no major mergers at z < 2, and with a low and constant gas accretion rate that 
keeps a stable angular-momentum direction. More violent merger or gas accretion histories lead to 
galaxies with more prominent bulges. Most disk galaxies have a bulge Sersic index n < 2. The galaxies 
with the highest bulge Sersic index tend to have histories of intense gas accretion and disk instability 
rather than active mergers. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: high-redshift — galaxies: 
spiral — galaxies: bulges — galaxies: interactions 



1. INTRODUCTION 

Spiral galaxy formation is shaped by a number of key 
dynamical processes, including galaxy-galaxy mergers, 
gas accretion from cosmic filaments and internal pro- 
cesses driven by disk instabilities. Understanding how 
these mechanisms combine to produce the range of mor- 
phologies observed in the Local Universe is one of the 
great challenges for galaxy formation theory. 

Mergers have been shown to be important for the for- 
mation of elliptical galaxies and for the growth of bulges 
in spirals. In particu lar, major mergers can transform 
disks into ellipticals |Toomre fc Toonire||1972| |Barnes j 



1992; [Naab et al.||1999| |JNaab & Bu rkert 2003| |CJox etaL 
2006), especially when they involve two relatively gas- 
poor galaxies, which can often be the case at low redshift. 
The end-product of gas-rich mergers is more complicated, 
with early studies arguing that the merger remnant can 



still contain a large rotating disk dSpringe l et al. 2005 



Robertson et al.||2 006; Hopkins et al. 2009 ) , while recent 



simulations resolving the star-forming regions and cap- 
turing supersoni c turbulence suggest that disk survival 
is very unlikely (Bournaud et al. 2011). 

While major mergers are violent events strongly im- 
pacting a galaxy's morphology, they are relatively rare 
in a ACDM context. Minor mergers (with mass ratios 
lower than 4:1) are mu ch more frequent and can also in- 
duce bulge formation ( Bournaud et al.| 2004| [2005 ; Cox 
et al.||2008[ ), ei ther by adding satell ite stars to the center 
of the galaxy ( |Aguerri et al.||200l| ) o r by triggering gas 
inflow s followed by central star bursts QEliche-Moral et ah] 
|2006| . However, the bulge fraction and Sersic index ol 
the remnant depend strongly on the detailed properties 
of the merger. Mass ratio is one of the key parameters: 
mergers with mass ratios between 4:1 and 1 0:1 are the 
most effic ient at heating the galactic disks (|Bournaud| 



et al. 2005). Smaller mergers only slightl y perturb the 
disks, except if they are frequent enough ( |Bournaud et al. 



2007bD 
IiT the end 



the effect of mergers on disk galaxies is 
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still debated, and even if they temporarily destroy disks, 
new disks could reform from the infall of fresh gas from 
cosmic filaments. At low redshift, mergers become less 
frequent and gas accretion can then play an important 
part in slowly building stellar disks around pre-existing 
bulges. 

Gas accretion however is not always slow and steady, 
particularly at high redshift where cold flows can rapidly 
accumulate large amounts of gas in galac tic disks and 
make them violently unstable ( jDekel et al. 2009b ). These 
gas-rich and unstable disks fragment and are observed as 
havin g a clumpy morphology in hig h-redshift galaxy sur- 
veys (El megreen et al.]|2004 2007). The clumps migrate 
toward the center of the galaxy under the effect of dy- 
namical friction (although it is debated if they survive 
the intens e star formation they host and the associated 
feedback - |Krumholz fc Dekel|2010| |Murray et al.|2010|) , 
where they participate in building the bulge, thus provid- 
i ng another channel for bulge formation at high redshift 



s 



Noguchi|199^|Immeli et al.|2004||Bournaud et al.|2007a[ 
jenzel et al.||20081 |Elmegreen et al.||20081 |Agertz et al. 



20Q9| |Dekel et al.||2009a| jCeverino et al ||2010 
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The picture becomes even more complex with the real- 
ization that local spiral galaxies host bulges with a whole 
range of morphologies. These are usually broken into 
two types, classical bulges and pseudo-bulges. Classical 
bulges have a high Sersic index (meaning their mass dis- 
tribution is highly concentrated) and a spheroidal shape. 
A mechanism explaining their formation is the violent 
relaxation of stars during mergers, either major or fre- 
quent eno ugh minor o nes (Bournaud et al.||2007b). Fur- 



thermore, Elmegreen et al.| ( |2008|) argue that clump coa- 
lescence in gas rich disks is actually similar to a merger in 
terms of orbit mixing and also pro duces classical bulges 
(although see |Inoue fc Saitoh||2012| showing it might not 
be the case). 

On the other hand, pseudo-bulges can be either flat- 
tened, with disk-like profile and kine matics or have a 
boxy-peanut shape (see the review by |Kormendy fc Ken-| 
|nicutt|[2004[ ). They are most likely the result of secular 
processes re-arranging disk material on timescales of a 
few billion years. Bars are in particular thought to play a 
central part in pseudo-bulge formation since they trigger 
gas inflows, resul ting in increased star formation in the 
centr al regions ([A thanassoula 2002} |2005 
20071 lEisher et al.| |2009| ). As bars grow 6. 



Heller et al. 

der 



vertical 

instabilities can give them the boxy /peanut shape ob- 
served in many edge-on spirals dCombes fc Sanders| l981 ; 
Combes et al.|1990||Raha et al.flOOll lE ureau fc Freeman 



1999; Debattista et al. 2006). Finally, minor mergers 
themselves are also a possible mechanism for triggering 
gas inflows and producing bulges with a low Sersic in dex 
flEliche-Moral et al.|[2006| |Mendez-Abreu et al.||2010 ). 

In the local Universe, pseudo- bulges are mostly found 
i n late- type galaxies, e specially in isolated environments 



(Durbala et al. 2008). The Milky- Way itself hosts a 



pseudo- bulge, with the contribution o f a classical bulg e 
limited to 8% of the total stellar mass (She n et al.|2 010). 
Even more stricking is the fact that 80% of local gal ax- 
ies (i.e. within the local 11 Mpc sphere) more massive 
than 1 9 are either bu lgeless or only have a pseudo 



bulge ([Fisher fc D rory 2011, see also |Kautsch et al.|2006 
and [Kormen dy et al ||2010p . The widespread existence 



of galaxies without a classical bulge is extremely puz- 
zling given the variety of mechanisms expected to pro- 
duce them, and the frequency of these mechanisms in our 
standard cosmological model. 

The formation of realistic spiral galaxies has indeed 
been a major challenge for cosmological simulations in 
the last 20 years. Early simulations suffered from ex- 
treme angular momentum loss during mergers, giving 
birth to galaxies with overly co ncentrated mass distri- 
butions and massive bulges (e.g., | Navarro fe Benz 



Navarro & White 1994 



Abadi et al. 1120030. T 



1991 



tie in- 



creased resolution ol recent simulations, as well as the 
inclusion of additional physics, in particular feedback 
from s upernovae, has considerably improved the situ- 



ation dGo vernato et al. 2007 Scannapieco et al. 2008 i 
|Pionte^fc Steinmetz 20 111. Zoom cosmological simu- 
lations can now produce Mi lky- Way mass galaxies with 
a bulge fraction of 0.2-0.3 (lAgertz et al. 2011 Guedes 
et al.| 2011 Brook et al.||2011a| , closer to that observed 



m our own Galaxy. In such simulations, the formation of 
galaxies with a small bulge relies on the use of a low star 
formation efficiency, a high threshold for star formation 
and/or gas recycling with galactic fountains. Although 
no clear consensus has been reached on how to model 
star formation and feedback, a picture emerges in which 
more disk-dominated galaxies form if their star formation 
histori es are delayed significant ly by the chosen subgrid 
recipe ( |Scannapieco et al.|2012| ). However, no simulation 
has ever been able to produce a Milky- Way mass bulge- 
less galaxy (note that this problem is partially solved 
for dwarf galaxies, where supernova feedback is much 
more efficient at removing baryons wit h a low angular 
momentum — see Governato et al.pOlO ). 

More generally, the intrinsic limitations of zoom simu- 
lations make it hard for them to study galaxy formation 
from a statistical point of view. In this type of sim- 
ulations, a given sub- volume of the box is simulated at 
high resolution, while the rest of the box is kept at a very 
coarse resolution. The high-resolution sub- volume is cen- 
tered on the galaxy that is being studied, and this volume 
needs to encompass all particles ending-up in the galaxy 
at z = 0. This can lead to a very complex, and possibly 
large, zoom volume, often making the technique highly 
inefficient, so that not all galaxies can be simulated this 
way. It is in particular difficult to simulate galaxies with 
high accretion and merger rates at all redshifts, although 
more isolated galaxies or massive galaxies with early as- 
semblies are more easily modeled. 

In any case, gathering a larger sample of such simula- 
tions is a challenge given today's supercomputers. Hence, 
it has been difficult to test the formation of galaxies of 
different Hubble types (at fixed physical recipes) , to link 
their morphology with their history or to study the con- 
nection between high reshift galaxy populations and their 
z = counterparts. 

In this paper we present a suite of 33 simulations per- 
formed with a n alternative zoom r e-simulation technique 
(presented in |Martig et al.|[2009[) . This technique con- 
sists of extracting the merger ana gas accretion histories 
of dark matter haloes in a large-scale cosmological sim- 
ulation and then re-simulating these histories at much 
higher resolution. Its main advantage is its low com- 
putational cost that makes it possible to gather a large 
sample of simulated galaxies and to follow their evolution 
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from z — 5 to z = 0. 

The simulated sample is made of 33 isolated galaxies 
with a halo mass between 2.7 x 10 11 and 2 x 10 12 M . 
At z = 0, we find galaxies with a large range of Hubble 
types, from bulgeless to bulge-dominated galaxies. Most 
of the galaxies host pseudo-bulges, with a Sersic index 
lower than 2. 

We study the evolution of their morphology with red- 
shift. We do not find any correlation between the mor- 
phology at z = 1 and at z = 0, with a whole range of 
possibilities for the z = 1 progenitors of spiral galax- 
ies (interacting galaxies, bulge-dominated systems, pure 
disks, unstable disks...). By contrast, there is a much 
better morphological correlation between z = 0.5 and 
0, with spiral arms and bars being mostly in place at 
z = 0.5. 

Focussing on the formation histories of galaxies with 
B/T< 0.3 (typically Sb and later types, and correspond- 
ing to 16 galaxies of our sample), we find that the most 
disk-dominated of these galaxies have both an extremely 
quiet merger and gas accretion history. By contrast, 
more violent merger or gas accretion histories give birth 
to galaxies with more prominent bulges. We find that 
the galaxies with the highest bulge Sersic index at z = 
are not those with mergers but those with intense gas ac- 
cretion at z = 1 and either early bar formation or other 
disk instabilities. 

In Section 2 we describe the simulation technique and 
discuss the selection of the sample. In Section 3 we 
present the techniques we used to perform morphological 
decompositions, while in Section 4 we discuss the general 
properties of the sample at z = 0. Section 5 is devoted 
to the high redshift progenitors of spiral galaxies, where 
we present the evolution of their mass, star formation 
rate, size and morphology with redshift. Finally, Section 
6 is focused on our most disk-dominated galaxies (with 
a bulge fraction lower than 0.3) and characterizes their 
merger and gas accretion histories. 

2. SIMULATIONS 

We study the evolution of 33 simulated galaxies from 
z = 5 to z = usin g the zoom-in techn ique described 
in Appendix A of |Martig et al] ( |2009| ). This tech- 
nique consists of extracting merger and accretion histo- 
ries (and geometry) for a given halo in a A-CDM cosmo- 
logical simulation, and then re-simulating these histories 
at much higher resolution, replacing each halo by a real- 
istic galaxy containing gas, stars and dark matter. 

In this Section, we will briefly recall the main details 
of the simulation technique, present the characteristics 
of the model galaxies, and explain how the sample of 33 
re-simulated galaxies was chosen. 

2.1. Simulation technique 

The collisionless cosmological simulation we use as a 
first step in this work has been performe d with the Adap- 
tive Mesh Refinement code RAMSES ( |Teyssier| |2002| ) , 
in a box of comoving length equal to 20 h _i Mpc. It 
contains 512 3 dark matter particles with each a mass of 
6.9 x 10 6 M . The cosmology is set to A-CDM with 
a m =0.3, ^ A =0.7, H =70 km s" 1 Mpc -1 and a 8 =0.9. 

The dark matter halos are detected i n the cosmological 
simulation using the HOP algorithm ( |Eisenstein fc Hut 



1998). Particles that do not belong to a gravitationally 
bound halo are also taken into account; we will refer to 
them as "diffuse" mass accretion when they are accreted 
by the studied halo. 

We extract the merger and diffuse accretion histories 
of a halo by tracking halos and particles crossing a fixed 
spherical boundary drawn around the target halo. We 
record the position, velocity and spin of each incoming 
satellite as well as the date of the interaction. We then 
perform a new simulation, in which we now replace each 
halo of the cosmological simulation by a galaxy made up 
of gas, star and dark matter particles, and each diffuse 
particle of the cosmological simulation with a blob of 
lower-mass, higher-resolution gas and dark matter par- 
ticles (see a deta iled explanation of how this is done in 
Appendix A of |Martig et al.| |2009). They start interact- 
ing with the main galaxy following the orbital and spin 
parameters given by the cosmological simulation. The 
re-simulation starts at z = 5 inside a 800 kpc-large zoom 
area, and follows the evolution of the main galaxy down 
to z = 0. 

The code we use for th e re-simulation is the P a rticle- 
Mesh code described in |Bournaud fc Combes| ( |2002| 
2003). The density is computed thanks to a Cloud-ln- 
Cell interpolation and the Poisson equation is then solved 
using Fast Fourier Transforms. Time integration is made 
using a leapfrog algorithm with a time step of 1.5 Myr. 
The spatial resolution (gravitational softening) is 150 pc 
and the mass resolution (particle mass) is 1.5 x 10 4 M 
for gas, 7.5 x 10 4 M for stars, and 3 x 10 5 M for dark 
matter. In Appendix A we present a resolution test, 
with the spatial and mass resolutions higher by a factor 
of two and six respectively. We find very little difference 
between the high-resolution and standard runs. 

Gas dynamics is modeled with a sticky-particle scheme 
(with a similar resolution of 150 pc) and star f orma- 



tion is comp uted with a Schmidt-Kennicutt law (Ken 
|1998[ ) with an exponent of 1.5 and an efficiency 
The threshold for star formation is set at 0.03 



nicutt 

M pc -3 (i.e. one atom per cubic centimeter), which 
corresponds to t he minimal dens ity for diffuse atomic 
cloud formation (Elmegreen 2002). We include the en- 



ergy feedback from supernovae explosions, and use a ki- 
netic feedback scheme: 20% of the energy of the super- 
novae is distributed to neighbouring gas particles (within 
a radius of 70 pc) in the form of a radial velocity kick. 
We also include the continuous gas return from high-, 
intermediate - and low-mass s tars following the scheme 
proposed by jJun g wiert et al. (2001) and used in |Martig| 
fc Bournaud| ( |2010| ). We indeed found stellar mass loss to 
be an important ingredient for the formation of realistic 
late- type galaxie s in cosmological simulations (Mar tig fc| 



Bournaud 2010). 



2.2. Model galaxies 

Each halo recorded in the dark matter-only cosmolog- 
ical simulation is replaced with a realistic galaxy, made 
of a disk (gas and stars), a stellar bulge and a dark mat- 
ter halo. The disks have a Toomre profile, the bulge is 
a Plummer sphere and the halo follows a Burkert profile 
flBurkert|[l995l |Salucci fc Burkert| [2000). The core radii 
ol our dark matter halos a re chosen according to the scal- 
ing relations proposed by Salucci fc Burkert] ( 2000 ) and 



4 



the halos are truncated at their virial radius. 

T he stellar mass of a g alaxy is set according to Table 
7 of Moster et al. (2010) as a function of its halo mass 
and of the redshift at which the galaxy is introduced 
in the simulation. All model galaxies are assumed to be 
disk galaxies without a bulge, except for massive galaxies 
(Mhaio > lO n M0) at z < 1, in which case the mass of 
the bulge is set to 20% of the stellar mass. 

The gas content in the disk also varies according to 
redshift and total mass: at z < 1, the gas fraction (with 
respect to total baryonic mass) is 0.3 for small galaxies 
(Mhaio < 1O 11 M ) and 0.15 for massive galaxies. At 
higher redshifts, the gas fraction is chosen independently 
of the total mass, and is set to 0.5 for 1 < z < 3 and 0.7 
for z > 3. No hot gas halo is included, which should be a 
sensible approximation given that the mass of the model 
galaxies is always lower than the criti cal mass for virial 
shocks to be stable (~ 10 12 M^, s ee |Birnboim fc Dekel 
2003| and |Dekel fc Birnboim|2006[ ) . 



2.3. Advantages and limitations of the simulation 
technique 

The technique we developed was designed to overcome 
some of the limitations of standard zoom simulations. 
Such a re-simulation tec h nique has al r eady been used by 
Kazantzidis et al |J [2008| ), |Read et aLlfl2008| ) and |Villalo- 
bos &; Helmil ( |2008| ) to study the stability of a stellar 
disk undergoing a series of collisionless minor mergers. 
A major difference with our work is that these simula- 
tions include no gas component (and consider mergers 
only above a given mass ratio) so that their scope is lim- 
ited and they cannot be used to study most mechanisms 
governing galaxy evolution. 

The main advantage of the re-simulation technique is 
the low computation time Q compared to zoom cosmo- 
logical simulations, so that running a large number of 
simulations is easier and statistical studies are possible. 
Another advantage is that it is possible to simulate galax- 
ies with all types of merger histories, in contrast to zoom 
simulations that need all the progenitors to be in the sub- 
volume treated at high resolution. Some merger histories 
are thus for now impossible to simulate with a standard 
zoom technique, in particular if they involve a merger 
at z ~ with a massive galaxy, in which case the sub- 
volume on which one has to zoom becomes very large. 
Finally, the fact that we decouple the evolution of our 
main galaxy from the expansion of the Universe keeps 
the physical resolution constant as a function of time at 
no additional cost (all positions are expressed in phys- 
ical units rather than in comoving units). This is not 
necessarily the case in cosmological simulations, that of- 
ten have a constant resolution in comoving coordinates, 
so that the physical resolution decreases with redshift 
unless (in the case of AMR simulations) additional levels 
of refinement are added, with a large computational cost. 

Of course there are also some major drawbacks. The 
most important is probably the large number of free pa- 
rameters for the model galaxies. Even if we carefully se- 
lect these parameter to be close to the observed galaxies, 
we will never reach the level of diversity found both in 

1 around 5 days for each simulation presented in this paper 




10 15 
x [Mpc] 




10 15 
x [Mpc] 

Fig. 1. — Location of the chosen galaxies within the cosmological 
box. The two panels show two different projections (in the xy and 
xz planes), each black dot corresponds to a dark matter halo (only 
halos more massive than 10 9 Mq are shown) and the red dots 
correspond to the 33 galaxies chosen for the re-simulation. The 
chosen galaxies sample the whole volume but avoid the densest 
regions. 

observations and in fully cosmological simulations. The 
case of our z = 5 galaxies is even more complicated, with 
very few data both from observations and previous sim- 
ulations. However, at z = 5 only a minor fraction of the 
final baryonic mass of a simulated galaxy is already in 
place: on average, we find that only 5% of the z = 2 
baryons are in place at z = 5. The mergers and insta- 
bilities happening between z = 5 and z = 2 redistribute 
these initial baryons, so that the initial assumptions do 
not make a big difference for our purpose here (see a 
test of a change of initial conditions in Appendix A of 
|Martig et aLl|2QQ9|> . We are thus careful to "discard" 
the early phase of the simulations and only analyse them 
after z = 2 to erase the memory of the initial conditions. 

Finally, an important limitation is linked with the 
sticky-particle model, that poorly treats the ho t gas 
phase (see also the discussion in Appendix B of |Mar-| 
|tig et aL]|2009| ). This technique is thus unable to treat 
the case of massive halos (above the critical mass of 
~ 10 12 M for virial shocks to be stable), for instance 
massive early-type galaxies, but also groups and clusters. 
We are limited to galaxies below the threshold for virial 
shocks, for which we know the cold accretion mode is 
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Fig. 2. — Examples of i-band surface brightness maps (face-on and edge-on, 70x70 kpc, color scale from 16 to 27 mag arcsec -2 ) for 2 
simulated galaxies, and the corresponding stellar mass growth histories (stellar mass computed within the optical radius). The two cases 
shown here represent two extremes: a disk-dominated galaxy with an extremely quiet history, and a bulge-dominated one with a much 
more violent history (two major mergers at z < 1). 

fected by th e relatively small box s ize of the cosmological 
simulation flPower fc Knebe||2006| ) . 

Appendix |B| shows lor each simulated galaxy several 
i-band surface brightness maps (computed using PE- 
GASE, with a Kroupa IMF, and without dust extinc- 
tion) at various redshifts, the i-band surface brightness 
profile at z = as well as the stellar mass evolution with 
time. The mass growth histories are usually very quiet 
for these galaxies, resulting in a disk+bulge (and possibly 
bar) morphology for all of them. Figure [2] shows exam- 
ples of i-band images and stellar mass evolution plots for 
two extreme cases: a disk-dominated galaxy (left panels) 
with an extremely quiet history, and a bulge-dominated 
one that has undergone several mergers at z < 1. 

We find that at z = 2 the simulated galaxies have at 
most reached 20% of their final stellar mass, at z = 1 
they have typically between 30 and 60 % of their final 
mass, and between 50 and 80% for z = 0.5. 

The redshift at which half of the final stellar mass is in 
place (zi/2) is shown in Figure [3] While Zi/ 2 is greater 
than 1 for 22% of the simulatea galaxies, the majority 
(69%) of galaxies reach half of their final stellar mass 
between z = 1 and z = 0.5. These results suggest a late 
assembly of our sample galaxies. We find one extreme 
case of a galaxy with only 0.001% of its mass in place 
at z = 2, and 4% at z = 1; this is a galaxy undergoing 
several major mergers after z = 1 and having a bulge- 
dominated morphology at z = (this is the galaxy shown 
on the right in Figure [2J. 

3. BULGE/DISK/BAR DECOMPOSITION 

We use two different methods to characterize the mor- 
phology of the simulated galaxies. One of them is based 
on stellar kinematics. It is usually robust but fails at dis- 
criminating bars and bul ges. The other techniqu e con- 



FlG. 3. — Redshift at which galaxies reach half of their final stellar 
mass ( £1/2 5 or formation redshift) as a function of their stellar mass 
at z = 0. The formation redshift of most galaxies is between 0.5 
and 1, suggesting a late assembly of their stellar content. 

predominant. In that case, the feeding of cold gas to 
the galaxy is most easily tied to the properties of the 
cosm ic-web filaments to which the galaxy is connected 
(e.g., |Dekel et al.||2009b| ). 



2.4. Sample selection 

We select 33 halos with a mass between 2.7 x 10 11 and 
2 x 10 12 M and that are relatively isolated at z = 0: no 
halo more massive than half of their mass can be found 
within 2 Mpc, and they are at least 6 Mpc away from 
one of the four most massive halos in the simulation box. 
The chosen halos are distributed across the whole volume 
of the simulation box but avoid the densest regions (see 
Figure [lj. We computed the local density contrast at z = 
in spheres of radius 8 Mpc around each halo, S = 
and find values ranging from -0.29 to -0.01, with all but 5 
galaxies in the range [-0.29,-0.22]. This confirms that the 
simulated galaxies are in underdense environments, with 
densities inter mediate between the Local Group and th e 
center of voids ( Romano-Diaz fc van de Weygaert|2 007). 

No additional criterion is set, in particular there is no 
constraint on the merger history of the chosen halos. The 
chosen halos represent 26% of the total number of halos 
in the simulation box within this mass range, and 69% of 
the halos satisfying both the mass and environment crite- 
ria. Note that in the mass range that we are considering 
the abundance and internal properties of halos are unaf- 



sists in using GALFIT (iPeng et al. 2002, 2010) to fit 



mock images of the galaxies with a sum ofHSulge, bar 
and disk profiles. This technique is less robust, especially 
for galaxies that are difficult to fit with smooth profiles, 
which is often the case at high redshift and which is the 
reason why we will only use this technique for the galax- 
ies at z = 0. On the other hand, GALFIT has the ad- 
vantage of giving measures that are easy to compare to 
observations. It also allows us to discriminate bars and 
bulges, and to measure their Sersic index. A detailed 
comparison of the relative advantages and limitations of 
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the tw o techniques can be found in [Scannapi eco et al. 
([2010). 



3.1. Kinematic decomposition 

We set the z-axis to the direction of the total an- 
gular momentum of all stars within the optical radius 
(R25 in g band). We compute the circularity parameter 
e = . Jz , m for each of the stellar particles within the op- 

Jcirc y-tt) 

tical radius, where j C irc(R) is the angular momentum a 
particle would have if it were on a circular orbit at its cur- 
rent radius. In a standard spiral galaxy, the distribution 
of e shows two peaks, one at ~ 1 corresponding to the 
rotating disk component, and one at e s ~ correspond- 
ing to the non-rotating spheroidal component (see Panel 
a in Figure [4| . We identify the spheroidal component by 
assuming it corresponds to a distribution of e symmet- 
ric around e s so that its mass is M s = 2 x M(e < e 3 ). 
Note that in some cases, particularly if the galaxy hosts 
a bar, e s is not 0, but most often between 0.2 and 0.4, 
which means the central regions of the galaxy have a net 
rotation movement (see Panel b in Figure [4]). 

In some cases however, e s can be difficult to measure 
if the two peaks are not clearly separated (see Panels c 
and d in Figure [4]). In this situation, we either set e s to 
the value corresponding to the peak of the distribution 
of e restricted to stars within the central inner 2 kpc of 
the galaxy (see Panel c), or if this does not help e s is set 
to (see Panel d). 

3.2. Photometric decomposition with GALFIT 

We create i-band images of the z = galaxies, extend- 
ing out to 1.5 times the optical radius. We mask out of 
the image satellite galaxies as well as any bright clump 
that would make the identification of the main compo- 
nents more difficult. We use GALFIT in two steps: first 
we try a bulge+disk model and then a bulge+bar+disk 
model (see some examples of decompositions in Figure 
[5]). Bulges and bars are assumed to follow a Sersic profile 
while disks are assumed to have an exponential profile. 
In a bulge+bar+disk decomposition, the bar is chosen 
to be the Sersic component with the lowest Sersic index, 
largest effective radius and smallest axis ratio. 

Between a bulge+disk and a bulge+bar+disk decom- 
position, x 2 i s often not sufficiently different to clearly 
identify the best model. In this case we add two ad- 
ditional criteria: the amplitude of the residuals both in 
the 2D image and in the radial profile. We then set the 
total bulge and/or bar luminosity to the corresponding 
Sersic fits. As far as the disk is concerned, taking the 
luminosity of the smooth exponential profile is not al- 
ways relevant, especially in the case of rings or lopsided 
features. We then measure the disk luminosity by sub- 
tracting the bulge and/or bar profiles to the total light 
distribution. 

Note that 5 galaxies present such a complex struc- 
ture that no satisfying GALFIT decomposition can be 
achieved, we remove them from studies involving photo- 
metric decompositions. 

3.3. Comparing the two techniques 

Figure [6] compares the disk-to-total (D/T) ratio ob- 
tained for each galaxy with the kinematic and photo- 
metric techniques. There does not seem to be any strong 



systematic offset between the values obtained by the two 
methods, although the scatter is quite large. For galaxies 
without a bar however, we find that GALFIT tends to 
give higher values of D/T compared with the k inematic 
decomposit i on. T his trend was also found by IScanna-| 
pieco et al. (2010). They argue that the most likely rea- 
son lor this difference is the fact that the photometric 
technique assumes that the disk exponential profile ex- 
tends all the way to the center of the galaxy, while the 
stars in the center most likely are dispersion-dominated 
and are instead attributed to the bulge component by the 
kinematic decomposition. Additional differences can be 
attributed to measuring the morphology from the mass 
distribution instead of the luminosity distribution. 

The fact that there is no systematic offset between the 
disk fraction obtained by the two methods suggests how- 
ever that the "spheroid" component identified kinemat- 
ically roughly corresponds to the bulge and bar compo- 
nents identified by GALFIT. 

4. PROPERTIES OF THE SAMPLE AT Z = 

We study the global properties of our simulated galax- 
ies at z = 0. We use in this section the photometric mea- 
surements to facilitate comparison s with observations, 
especially with Weinzirl et al. (2009), who studied a sam- 
ple of local massive spirals with a GALFIT decomposi- 
tion technique similar to the one adopted here. We show 
in Figure [7] the cumulative distribution of the bulge-to- 
total ratio as measured with GALFIT 

The simulated galaxies span a large range of bulge frac- 
tions, from nearly bulgeless galaxies with B/T of the or- 
der of 0.05 to bulge-dominated galaxies with B/T close 
to 0.8. However, the distribution of the bulge fractions 
differs significantly from the observed one: B/T tends 
to be too high in the simulations. Figure [7] shows that 
only 32% of our simulated galaxies hav e B/T< 0.2, a 
value much lower than the 66% found by Weinzirl et al. 
(2009). Similarly, 40% of the observed sample is found 
at B/T< 0.1, a value that is only reached by 14% of 
our sample (these nearly bulgeless galaxies are shown in 
Figure [25]). In spite of the success of the simulations 
at producing galaxies with small bulges, there is thus 
strong evidence that not enough simulated galaxies are 
produced with a low B /T. Although the images used by 
Weinzirl et al. (2009) are in H-band while ours are in 
i- band and we do not include dust extinction, and while 
these differences will affec t the measured bulge fractions 
(Graham & Worley 2008), this is probably not enough 
to reconcile the simulations with the observations. 

In addition, even if we probe similar stellar mass 
ranges, our simulated galaxies are likely more isolated 
than the observed sample, that is only magnitude- 
limited, with no restriction on environment. However, 
if we compare the simulations to samples of isolated 
galaxies, the discrepancy becomes even greater: the ob- 
served population of isolat ed galaxies is strongly domi- 
nated by la te-type spirals ( |Sulentic et al.|[2006| |DurbaIa 
eT^|2008l ). 

Note also that we do not find any pure elliptical galax- 
ies, but this is not so worrying since they ar e rather 
rare, especially in isolated environments (e.g 
19801 iHogg et al.||2003| |Croton et aL]|20Q5|). 
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&; White| ( |201l| ) perform bulge/ disk decompositions ol 
SDSS galaxies and also show that very few galaxies ac- 
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Fig. 4. — Examples of bulge/disk decompositions based on kinematics. Each panel illustrates a different possible case; for each galaxy 
the distribution of e = j z /jdrc is shown (together with the resulting spheroid/disk decomposition), as well as i-band face-on and edge-on 
images (50 kpc x 50 kpc images). Panels a and b correspond to cases where two peaks are easily identified in the distribution of e so that 
e s is straightforward to measure (note that Panel b shows a galaxy with a bar, for which e s ~ 0.3). Panel c and d correspond to galaxies 
for which the spheroidal component is harder to identify. In this situation, we either set e s to the value corresponding to the peak of the 
distribution of e restricted to stars within the central inner 2 kpc of the galaxy (Panel c), or if this does not help e s is set to (Panel d). 
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Fig. 5. — Examples of GALFIT decompositions. The top panels 
show a bulge/disk decomposition (left: the simulated galaxy, mid- 
dle: the GALFIT model, right: the residuals), while the bottom 
panels correspond to a bulge/bar/disk decomposition. 

tually have bulge fractions greater than 0.8. 

Another difference with observations is that we do not 
find any correlation between stellar mass and bulge con- 
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tent (Figure pi), contrary for instance to 
(2009) that si low a trend to higher bulge 
higher stellar masses. This difference might be insignifi- 
cant due to our small sample and due to different meth- 
ods used to measure stellar masses in observations and 



(J2009) estimate the stellar 



simulations ( Weinzirl et al. 

mass of a galaxy from its B — V color, they estimate typ- 
ical errors to be within a factor of 2-3), but it could also 
reflect a more profound discrepancy between our simula- 
tions and observations. 



Fig. 6. — Comparison of photometric and kinematic decompo- 
sitions. We show the disk-to-total ratio (D/T) obtained with the 
kinematics as a function of that obtained with GALFIT. Barred 
galaxies correspond to the red dots (we define barred galaxies as 
having a bar-to-total ratio greater than 0.1 according to GALFIT), 
non-barred galaxies to the blue dots. 

Figure [9] shows the bulge Sersic index as a function of 
the bulge Traction. Most of our simulated galaxies have a 
low Sersic index, the majority lower than 2, indicating a 
large fraction of pseudo-bulges within our sample. This 
is consist ent with the obse r ved p roperties of local disk 
galaxies; Laurikainen et al. (2007) find for instance that 
index of bulges is on 



the Sersic index of bulges is on a verage lower than 2 
for all morphologica l types (see also |Durbala et al.| 2008 
Weinzirl et aT|2009 ). 
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Fig. 7. — Cumulative distribution of the bulge-to- total ratio mea- 
sured with GALFIT, compared to observations by Weinzirl et al. 
(2009). The simulations explore a large range of bulge fractions, 
from nearly bulgeless to bulge-dominated galaxies. However, the 
fraction of galaxies with B/T< 0.2 is only 0.34 for our sample of 
simulations, vs a f raction of 0.66 for th e sample of local massive 
spirals observed by Weinzirl et al. (2009). This is an important dis- 
crepancy, even if our sample is small and statistical comparisons 
are hard to perform. 
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Fig. 8. — Bulge-to-total ratio (measured with GALFIT) as 
a function of stellar mass at z = for the simulated galax- 
ies. We do not find any trend between bulge content and stel- 
lar mass, contrary to observations showing that more massive 
galaxies tend to be more bulge-dominated (e.g., | Weinzirl et al.| 
|2009] ). 



The bar content of our galaxies also appears roughly 
consistent with observations: 70% of our galaxies host a 
bar, to be compared to values of 60- 70% for local galaxies 
when observed in the nea r- infrared (Eskridge et al. 2000; 
jMarinova fc Jogee 2007). The simulated galaxies span 



a large range ol bar tractions, from galaxies with no bar 
to galaxies where the bar luminosity amounts to nearly 
half of the total luminosity. In a companion paper, we 
examine the properties of bars in our sample, including 
bar formation and the redshift evolution of the fraction 
of barred galaxies. 

To summarize our results thus far: 

• Simulated Milky Way- mass spiral galaxies show at 
z = a wide range of bulge fractions, from bulge- 
less to bulge-dominated galaxies, though the bul- 
geless disks are under-represented. 

• These galaxies also display a diverse range of bar 
fractions 
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Fig. 9. — Bulge Sersic index as a function of the bulge-to-total 
ratio measured with GALFIT. Most of the simulated galaxies have 
a bulge Sersic index lower than 2, suggesting a high incidence of 
pseudo-bulges, which is consist ent with the observed properties of 
local spiral s (iLaurikainen et al.|2007||Durbala et al. 2008 Weinzirll 

|et al.|2 009). — 

* The majority of bulges in our sample would be 
classified as pseudo-bulges and not classical bulges 
based on their Sersic index. 

Our simulations fail to reproduce the observed distri- 
bution of bulge fractions (with too few bulgeless galaxies 
within the simulated sample), as well as the observed 
correlation between bulge fraction and stellar mass. In 
spite of these failures, the sample of simulated galax- 
ies remains relevant to study the redshift evolution of 
galaxy morphologies and the physical mechanisms build- 
ing bulges. The next sections will investigate the nature 
of the high-redshift progenitors of the simulated galaxies, 
and will also focus on the most disk-dominated cases at 
z = and study their formation histories. 

5. THE HIGH-REDSHIFT PROGENITORS OF SPIRAL 
GALAXIES 

In this section, we study the evolution with redshift 
of the star formation rates, gas fractions, sizes and mor- 
phologies of the simulated galaxies. Figure [10] shows i- 
band surface brightness maps for a subset oTsimulated 
galaxies at z = 2, 1, 0.5 and 0. The maps for the whole 
sample can be found in Appendix [Bj 

5.1. Mass growth and star formation histories 



Figure 11 shows the star formation rate of the simu- 
lated galaxies as a function of their stellar mass at red- 
shift 2, 1 and 0. The simulation s are compared to obser- 
vations of s tar-forming galaxies ( Elbaz et al.|2007 



Daddi 



eTaLl[2007| ). We find an agreement between simulations 
and observations at z = 2 and 0, while at z = 1 the 
simulated galaxies have systematically lower SFRs than 
the observed ones. This disagreement could be (at least 
partially) due to selection effects, for instance it could 
be that z = spirals were not on the star-forming se- 
quence at z — 1, i.e. that the galaxies on the star form- 
ing sequence at z = 1 have evolved into bulge-dominated 
galaxies at z = 0. However, another explanation is a 
too short gas consumption timescale in the simulation 
between z = 1 and 2 (maybe linked with an overly weak 
feedback implementation), which would result in low gas 
fractions at z = 1, and thus low SFRs. 

To investigate this issue, Figure [12] shows the evolution 
with redshift of the gas fraction, denned as the total gas 
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Fig. 10. — Morphological evolution from z = 2 to for a subset 
of simulated galaxies. We show i-band surface brightness maps 
(face-on projection, 70x70 kpc) at z= 2, 1, 0.5 and 0. 

mass within R25 divided by the total mass of gas and 
stars within R25. While the values for z = 2 seem in 
agreement with observations, the z = 1 gas fractions 
are indeed probably a bit too low (|Daddi et al. 2010; 
Tacconi et al.|[20T0 Geach et al.||201l| ), especially since 
they reflect the whole gas content of the galaxies and not 
only their molecular gas content (precise comparisons are 
thus difficult, in addition to observational uncertainties 
in the CO-to-H2 conversion factor). It is interesting to 
notice a tight sequence of decreasing gas fraction both 
with decreasing redshift and increasing stellar mass. 



5.2. Sizes 

We study the half-mass radius of the galaxies as a func- 
tion of their stellar mass. Figure [l3| shows the results at 
z = 2 , 1 an d 0, comp ared with observa tions by Shen 



aD (|2003| ) for z = 0. |Shen et al.|j2003| ) actually com- 
ite K50 in the r-band, which is different from using the 



et 

pute 

mass distribution. We find however that our simulations 
agree nicely with their measures. The simulated galaxies 
have reasonable sizes, a fact that had been an issue for 
galaxy formation simulations in the last few decades but 
not so much now (|Brooks et al.|2011 |Guedes et al.|2011 ; 
Brook et al. 2011 bp . We note that at lixed stellar mass 
our disk- dominated galaxies tend to have larger values of 
R50 with respect to the bulge-dominated ones, which is 
also found in z = observations. 

A difficulty when trying to compare with observations 
at higher redshifts comes from the variety of morpholog- 
ical types and the limited mass range that we explore. 
For instance studies of the sizes in the zCOSMOS sam- 
ple at z ~ 1 ( |Sargent et al.|2007||Maier et al.|2QQ9| > focus 
on galaxies with a stellar mass greater than 5 x 1O 1O M , 
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Fig. 11. — Star formation rate as a function of stellar mass for 
the simulated galaxies at z = 2 (blue dots), 1 (red triangles) and 
(green squares), compa red to the correla tions observed f or star 
forming ga laxies (z = 2, Daddi et al. 2007 blue line; z = 1, Elbaz 
et al. 2007 red line; z = U, Elbaz et al. 2UU7 green line). The rilled 
symbols correspond to galaxies that are disk-dominated at z = 
(i.e. with a photometric B/T lower than 0.3), the empty symbols 
correspond to the rest of the sample. 
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Fig. 12. — Gas fraction as a function of stellar mass for the sim- 
ulated galaxies at z = 2, 1 and 0. The gas fraction is denned as 
the total gas mass within R25 divided by the total baryonic mass 
within R25. The filled symbols correspond to galaxies that are 
disk-dominated at z = (i.e. with a photometric B/T lower than 
0.3), the empty symbols correspond to the rest of the sample. 
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Fig. 13. — Half-mass radius as a function of stellar mass for the 
simulated galaxies at z = 2, 1 and 0. The filled symbols correspond 
to galaxies that are disk-dominated at z = (i.e. with a photo- 
metric B/T lower than 0.3), the empty symbols correspond to the 
rest of the sample. The z = o bserved relations for late-type and 
early- type galaxies ( |Shen et al.|2003) are also plotted. 
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a threshold that is only reached by 3 of our simulated at 
this redshift. O ur results are thus no t necessarily in con- 
tradiction with Sargent et al. (2007), who find as many 
galaxies with a radius between 5 and 7 kpc at z = 1 
and z = 0, whereas we find many more of these galaxies 
at z = than at z = 1. Additional discrepancies could 
come from using light profiles instead of mass profiles 
(and dust extinction could play a part: if extinction is 
underestimated in the central regions, the half-light radii 
could be overestimated). 

Overall, we find that the simulated galaxies at z = 2 
and 1 seem to follow t he same relation as the local galax- 
ies ( |Shen et al. 2003 for the SDSS): disk galaxies grow 

This idea 



along a sequence m the radi us-mass plane . 
is supporte d by observations (Bar den et al.||2QQ5 Maier 



et al. 2009), and it has alrea dy been seen in other sim- 
ulations (Bro oks et al.| |20TTj) and semi-analytic models 
Somerville~et al. 1 120081 iFirr 
utton et al.||2010 
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IN one ol our high redshift galaxies are particularly com- 
pact, which does not allow us to draw any conclusion 
on the nature and frequency of compact ellipticals. This 
suggests however that these galaxies are not found within 
the progenitors of isolated, Milky Way-mass galaxies 
(early- type galaxies are found in denser environments), 
and that instead they follow a different evolutionary 
track. 

5.3. Morphologies 

To characterize the morphological evolution of the sim- 
ulated galaxies, we first study bulge formation as a func- 
tion of redshift. Figure [l4| shows, as a function of the final 
bulge mass, the fraction of the bulge already in place at 
z = 2, 1 and 0.5. We use the kinematic decomposition 
so that "bulge" here actually means "bulge+bar" (as ex- 
plained in Section 3, using GALFIT is extremely difficult 
on high redshift galaxies). The color coding in Figure 14 
indicates the value of e s at z = for each galaxy ancTTs 
related to the amount of rotation in the central regions, 
i.e. to the presen ce o f a bar or a rotating bulge (e s is 
defined in Section [XT] as the average circularity for stars 
in the central spheroidal component, it is equal to for 
a non-rotating component). 

We find that only a minor fraction of the bulges were 
already in place at z = 2: galaxies typically have at most 
15% of their bulge-bar formed at this redshift. A strong 
evolution occurs between z = 2 and z = 1, and our 
simulations show a few cases of galaxies that have their 
entire z = bulge in place by z = 1; these tend to be the 
galaxies without bar at z = (blue colors on Figure [l4| ) . 
Galaxies with a strong bar at z = have by contrast still 
a lot of bulge-bar growth at z < 0.5. Note that Figure [l4| 
shows a few cases of galaxies with a greater bulge mass 
at z = 0.5 or 1 than at z = 0. This can be due either 
to errors in bulge mass determinati on or to a transfer of 
mass from the bulge to the disk (see Martig & Bournaud 
|2010| for an example of how stellar mass loss from bulge 
stars can transfer mass from a bulge to a disk). 

Instead of the total bulge mass, Figures [15] and [l6| show 
how the bulge fraction of the simulated galaxies evolves 
with redshift. Figure 15 studies the (potential) correla- 
tions between B/T at z = 2, 1, 0.5 and for the whole 
sample, while Figure 16 shows some examples of redshift 
evolution of B/T for a subset of galaxies. We find that 



at z = 2 (left panel of Figure 15), most galaxies (except 
one) are bulge-dominated according to their kinematics 
(they all have a relatively low mass, however, so that 
the measures are noisy). Many changes happen between 
z — 2 and z = 1, and by z = 1, the whole range of B/T 
is explored, from pure bulges to pure disks, with a few 
cases of mergers and of disks undergoing a phase of vio- 
lent instabilities (see the i-band surface brightness maps 
in Figure 10 and in Appendix [B|. Interestingly enough, 
no correlation is found between the morphology at z = 1 
and 0: pure bulges and pure disks at z = 1 can end up 
with the same bulge content at z = 0. This is also il- 
lustrated by Figure 16 that shows the evolution of B/T 



for ten galaxies endmg-up with a similar bulge fraction 
at z = 0, and that shows the large diversity of possible 
evolutionary tracks. 

Finally, we divide the z = 1 galaxies into bulge- 
dominated and disk-dominated samples, using the me- 
dian B/T at that redshift (equal to 0.54) as the limit 
between the samples. For each sample, we show in Fig- 
ure fl7| the range of final values of B/T reached at z = 0. 
Botnoulge- and disk-dominated z = 1 galaxies have a 
large (and similar) range of z = morphologies. We 
might even see a slight anti-correlation between the bulge 
content at z = 1 and 0, with z = 1 disk-dominated galax- 
ies tending to be more bulge-dominated at z = 0. 

This is also related to the fact that pure disks at z = 1 
do not remain bulgeless until z = 0: either a bar or 
a bulge grows in between. These pure disks are indeed 
very unstable, and even a small perturbation can destroy 
them. We find for instance the case of a galaxy that 
is bulgeless until z = 0.2, when a fly- by happens and 
trig gers bulge formation (this is the last galaxy of Figure 
|32| ). This is c onsistent w i th the studies of minor mergers 
performed by |Cox et al.| ( |2008| , where bulges are shown 
to stabilize disks, and to suppress merger-driven inflows 
and associated star formation, so that bulgeless galaxies 
are much more fragile. 

In contrast to the high redshift results, there is more 
correlation between B/T at z = 0.5 and z = (right 
panel of Figure 15). This correlation is actually closer 
to a bimodal behaviour since galaxies appear to be sepa- 
rated into two different sequences. For one sequence, the 
bulge-bar content decreases between z = 0.5 and 0, these 
are the galaxies that have no (or nearly no) bar at z = 
(blue circles in Figure 15). The galaxies on the other se- 
quence show an increase of their bulge-bar content, and 
these are galaxies having a substantial bar at z = (yel- 
low to red circles). This suggest a scenario for which, in 
galaxies without a bar, the bulge is in place by z = 1-0.5, 
and some important disk growth happens at z < 0.5. By 
contrast, in galaxies that develop a bar, there is some 
bar growth between z = 0.5 and 0, accompanied with 
late bulge formation. 

A visual inspection of the face-on images in Appendix 
[B] shows that bars and spiral arms are present in a few 
galaxies at z = 1, but are rare. They become much more 
common at z = 0.5. 

To summarize this Section: 

• simulated galaxies tend to follow the observed scal- 
ing relations between stellar mass, SFR, gas frac- 
tion and sizes. One exception could be a too strong 
gas consumption between z = 2 and 1, leading to 
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Fig. 14. — Evolution of the mass of the dispersion-dominated component (bulge+bar as measured by kinematics technique) with redshift. 
The three panels show the fraction of the mass of this component already in place at high redshift (left panel: z = 2; middle: 2 = 1; right: 
z = 0.5) as a function of the total bulge+bar mass at z = 0. The color code indicates the value of e s at z = for each galaxy, it is thus 
related to the amount of rotation in the central regions, i.e. to the presence of a bar or a rotating bulge. 
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Fig. 15. — Evolution of the (bulge+bar)-to-total ratio (bulge+bar as measured by kinematics) with redshift. The three panels show B/T 
at different redshifts (left panel: z = 2; middle: 2 = 1; right: z = 0.5) as a function of B/T at z = 0. The color code is the same as in 
Figure [14] 

too low SFRs and gas fractions at z = 1 

• there is no correlation between the morphology at 
z = 1 and at z = 0, and there is a whole range 
of possibilities for the z = 1 progenitors of spirals 
galaxies 

• the morphology at z = 0.5 is much closer to the 
final morphology, with spiral arms and bars being 
mostly in place at z = 0.5 

• the main epoch for bulge formation seems to be 
around z = 1, and we find some late bulge growth 
that accompanies bar formation at z < 0.5 




6. FORMATION HISTORIES OF THE MOST 
DISK-DOMINATED GALAXIES 

In this Section, we discuss the merger and accretion 
histories leading to the formation of the most disk- 
dominated galaxies of our sample. The reason to focus on 
these galaxies is that they are usually rare in cosmologi- 
cal simulations, and we wish to understand what is spe- 
cial about them. We use GALFIT photometric decom- 
position throughout this Section to classify morphology 
because it is the only way we can discriminate between 
bulge and bar. We select galaxies with B/T< 0.3 at 
z = 0, which is usually used as a threshold for identifying 
galaxies of Hubble types later than Sb, and corresponds 
to 16 galaxies in our sample. 



Fig. 16. — Evolution of the mass of the (bulge+bar)-to-total ratio 
with redshift for galaxies with B/T< 0.4 at z = 0. A large diver- 
sity of histories can give birth to galaxies with a similar bulge-bar 
fraction at z = 0. The two red lines correspond to galaxies with 
the same final B/T, but at z = 1 one is a pure elliptical, the other 
is a pure disk. The general trend between z = 1 and z = is a 
decrease in the bulge fraction, except for galaxies that are nearly 
bulge less at z = 1. 

6.1. Method 

We fix a boundary, which is a sphere centered on the 
main galaxy, and with a radius equal to 1.5 x R2$(z = 0). 
We follow the inflow of gas and stars through this bound- 
ary with a time resolution of 37.5 Myr. We compute both 
the mass and the angular momentum of the inflowing ma- 
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Fig. 17. — Distribution of the final values of the (bulge+bar)- 
to-total ratio for galaxies having at z = 1 either a low or a high 
B/T (the distinction between low and high B/T is made with re- 
spect to the median B/T at z = 1). Both bulge-dominated and 
disk-dominated z = 1 galaxies have a large range of z = descen- 
dants. 

terial. We discriminate between particles belonging to a 
satellite (that then correspond to mergers or fly-by s) and 
diffuse accretion of gas. 

We also compute this inflow only for particles that end- 
up within the optical radius at z = 0. This allows to 
discriminate fly-by s (where a satellite passes through the 
boundary but only leaves an insignificant amount of stars 
and gas in the main galaxy) and mergers (where a sig- 
nificant amount of mass from the satellite ends-up in the 
final galaxy at z = 0). 

6.2. Merger histories 

To study the merger histories of the simulated galaxies, 
we will make a distinction between major (mass ratios 
from 1:1 to 1:4), intermediate (1:4 to 1:10) and minor 
(smaller t han 1:10) merge r s. Th i s dist inction follows the 



work by Bournaud et al. (2004 2005) showing that in- 



termediate mergers can significantly affect galaxies, pro- 
ducing remnants with spiral-like luminosity profiles but 
with elliptical-like kinematics. 

Figure [18] shows the contribution of major and inter- 
mediate mergers to the baryonic growth of the simulated 
galaxies (more exactly the fraction of gas and stars within 
the optical radius at z = that have been brought in by 
these mergers) as a function of their final bulge fraction. 
The top panel of this Figure takes into account mergers 
at all redshifts. We find an overall low contribution (no 
more than 40%) of major and intermediate mergers to 
the mass growth of these galaxies, and a slight trend of 
an increasing importance of mergers with increasing final 
B/T. The galaxies with the lowest z = bulge fractions 
all assemble only 10-20% of their mass through major 
and intermediate mergers. Overall, the values we find 
are in rough agreement with other cosmological simula- 
tions studying the contrib utions of mergers and accretion 
to the growth of galaxies jMurali et al . 2002 ; Semehn~fc| 
|Combes|[2005l |Dekel et al.|2009b[ |L ? Hu'llier et al.||2012p . 
The bottom panel oi Figure shows only the contribu- 
tion of mergers occurring at z < 2 (this is when most of 
the mass growth happens). A large majority of galaxies 
assemble less than 20% of their baryons through major 
and intermediate mergers at z < 2, and all but one of 
the galaxies with B/T< 0.15 undergo only minor merg- 
ers after z = 2. 
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Fig. 18. — Fraction of mass (gas and stars) that joined the galax- 
ies through intermediate and major mergers as a function of their 
B/T at z = 0. We either consider mergers happening at all red- 
shifts (top panel) or only at z < 2 (bottom panel). 



Figure 19 shows the mass ratio of the largest merger 
undergoneby the galaxies after a given redshift (2, 1 and 
0.5) as a function of their final B/T. We find that some 
of the most bulge-dominated spirals of our sample have 
undergone at least one major merger between z = 2 and 
z = 1, but no galaxy has had a major merger at z < 1. 
Only 2 of the 16 galaxies have at least one intermediate 
merger at z < 1, and none of them has such a merger at 
z < 0.5. Note also that 7 galaxies (nearly half the cases) 
only experience minor mergers at z < 2. 

Focusing now on the left panel of this Figure, i.e., on 
the mass ratio of the largest merger after z = 2, one 
stricking feature is the large range of possible final bulge 
fractions for a similar merger history. Galaxies under- 
going only minor mergers at z < 2 can either be nearly 
bulgeless or have a bulge fraction of 0.2. 

The Sersic index of the bulges at z — (shown by the 
color code in Figure [l9|) can help to understand part of 
this messy situation. Three galaxies have a final bulge 
with a Sersic index greater than 2, and they appear 
located in the same region on the left panel of Figure 
19). These three galaxies only undergo minor mergers at 
z < 2, and have final bulge fractions of 0.10-0.2. They 
appear as a distinct population to galaxies who have also 
grown through minor mergers but instead have a lower 
B/T by z = 0. They are also distinct from bulge-rich 
galaxies with a more violent merger history. This actu- 
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Fig. 19. — Mass ratio of the largest merger undergone by the galaxies at z < 2 (left panel), 1 (middle) and 0.5 (right) as a function of 
their B/T at z = 0. The arrows mark lower limits. The color code indicates the Sersic index of the bulge at z = 0. 



ally suggests that mergers are not responsible for most 
of the bulge growth in these three galaxies, and we will 
get back to this in more detail in the next Section. Ex- 
cluding these three galaxies, and only keeping in mind 
those with a Sersic index lower than 2, we find a clearer 
trend of increasing bulge fraction and increasing Sersic 
index with increasing mass ratio of the largest merger 
undergone after z = 2. A large scatter still remains how- 
ever, and we do not find a one-to-one correlation between 
merger history and z — properties. 

An additional fact to notice here is that the galaxies 
undergoing a major merger at high redshift are found to 
have final Sersic indices in the intermediate range, with 
values between 0.6 and 1.5. These are relatively low val- 
ues for galaxies undergoing a major merger, and they do 
not fit in the standard picture of major mergers building 
classical bulges. We have studied the properties of these 
galaxies after the mergers (waiting for the remnant to be 
relaxed), and find that the bulges already have low Sersic 
indices: the low z = values are not the result of bulge 
evolution following (for instance) disk re-growth, rather 
the structure of the bulge is already in place after the 
merger. This might be due to the gas-rich nature of the 
mergers happening at high redshift. 

An interesting case is a galaxy undergoing a merger 
with a mass ratio 1:4.1 just after z = 1, but that ends 
up as extremely disk-dominated at z = (B/T=0.07, 
Bar/T=0.06). We find that the interaction increases the 
radius of the disk and triggers the formation of a spiral 
structure. This galaxy at z = has a very small bulge 
with a higher Sersic index (n=1.3) than the other B/T< 
0.1 galaxies (for which n is between 0.2 and 0.6), and a 
lower bar fraction. 

To summarize, we find that galaxies with B/T< 0.3 
have extremely quiet merger histories, both in terms of 
the mass ratio of most important merger they undergo 
and in terms of the fraction of baryons brought in by 
major and intermediate mergers. None of them have a 
major merger at z < 1, and only two of them have an 
intermediate merger after z = 1 (and these take place 
between z = 1 and 0.5). The most disk-dominated cases 
(B/T< 0.1) tend to have even quieter histories. If we 
exclude the three galaxies having a z = bulge with 
a Sersic index greater than 2, we find a correlation (al- 
though a very weak one) between the bulge fraction at 
z = and the mass ratio of the largest merger under- 
gone after z = 2. This suggests that mergers cannot 
be solely described by their mass ratio, but are highly 



complex events (mergers of a similar mass ratio can have 
very different consequences as a function of the orbit, gas 
fraction...). The existence of galaxies with quiet merger 
histories and large bulge fractions indicates that the mor- 
phology is partly determined by the gas accretion history, 
as we discuss in the next Section. 

6.3. Gas accretion histories 

We study the evolution with time of the average gas ac- 
cretion rate (Figure 20 ) and angular momentum (Figure 



21), and look for potential correlations with their bulge 
content at z = 0. Note that we only consider here the 
accretion of diffuse gas, and do not take into account gas 
brought in by satellite galaxies. 

We show in Figure [20] the average gas accretion rate 
for our simulated galaxies, separated in accretion at high 
(z > 1 — left panel) and low (z < 1 — middle panel) 
redshift s. The choice of z = 1 as a threshold is motivated 
by the results of Section |5.3| where we show that this is 
the typical epoch of bulge formation. When studying 
the influence of the gas accretion rate at z > 1 on the 
bulge fraction, we find that galaxies can be divided into 
two different populations: a tight sequence of increasing 
bulge fraction with increasing gas accretion rate, and a 
clearly distinct cloud of galaxies with low gas accretion 
rates but high bulge fractions (marked with blue and red 



points respectively in Figures 20 and 21). This cloud is 



actually made of galaxies that undergo at least a merger 
with a mass ratio greater than 1:5 after z = 2. If mergers 
are responsible for most of their bulge formation, it is 
thus not surprising to find them as outliers on this Figure. 
These galaxies also have a low gas accretion rate at z < 1 
(middle panel of Figure 20), so that disk re-building is 
harder, thus increasing tne chance for these galaxies to 
have a large bulge fraction at z = 0. 

Setting apart these 6 galaxies, we find a good correla- 
tion between the bulge content of a galaxy and its gas 
accretion history. The most disk-dominated cases have 
relatively low gas accretion rates, both at high and low 
redshift. The right panel in Figure 20 studies gas ac- 
cretion between z = 2 and 1, and shows the fraction of 
this accretion that happens in bursts. We define a burst 
as a period of gas accretion at a rate that would double 
(or more) the baryonic mass of the galaxy in less than 
1 Gyr. We find that the most disk dominated galaxies 
do not have any bursts of gas accretion between z = 1 
and 2, whereas galaxies that do have some bursts end up 
more bulge-dominated by z = 0. Note that these bursts 
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Fig. 20. — Gas accretion rates and their relation with the bulge fraction of the simulated galaxies. The left panel shows the average gas 
accretion rate at high redshift (z > 1), while the middle panel shows it for z < 1. The right panel shows for 1 < z < 2 the fraction of diffuse 
gas that is accreted in bursts, i.e. with a rate that would double (or more) the current baryonic mass of the galaxy within 1 Gyr. The most 
disk-dominated galaxies do not undergo any burst of gas accretion at these redshifts. The red dots highlight the cloud of galaxies with a 
low gas accretion rate at z > 1 but a high bulge fraction, that appear as outliers from the tight sequence found for the other galaxies. All 
these outliers are actually galaxies undergoing at least one merger with a mass ratio greater than 1:5 after z = 2, so that most of their 
bulge content is probably not determined by their gas accretion history. 
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Fig. 21. — Evolution with time of the direction of the angular momentum of the accreted gas. We study the evolution with time of cos(</>), 
where <j> is the angle between the angular momentum of the gas accreted at a given time and the final angular momentum of the galaxy. 
The left panel shows the mass- weighted average value of cos(</>) for z > 1, the middle panels shows the mass- weighted average cos(< />) f or 
z < 1, and the right panel shows the difference between these two values. The red dots correspond to the same galaxies as in Figure |20] 

of gas accretion are not necessarily related to mergers: 
Figure [20] shows cases of galaxies with a quiet merger 
history but a bursty gas accretion (the blue points with 
B/T> 0.15 on the right panel). 

Finally, we measure the angular momentum of the ac- 
creted gas, and how the direction of this angular mo- 
mentum changes with time. We define <p as the angle 
between the angular momentum of the gas accreted at a 
given time (at a constant radius of 1.5 x R2$(z = 0)) and 
the final angular momentum of the stellar component. 
We show in Figure 21 mass- weighted average values of 
cos((/>) at low and high redshift. We find that the galax- 
ies with the highest bulge fraction tend to have a smaller 
cos((/>) at high redshift, corresponding to a larger offset in 
the direction of the angular momentum of the accreted 
gas. At z < 1 most galaxies accrete gas with an angular 
momentum more closely aligned with the final z = mo- 
mentum (which is not surprising since this gas is often 
that which gives birth to a large fraction of the galactic 
stellar disk). The right panel of Figure 21 shows the dif- 
ference between cos((/>) at z > 1 and at z < 1. We find 
that the galaxies with the lowest bulge fraction undergo 
very little changes in cos((/>) as a function of time, and 
that large changes are linked with higher bulge fractions 
(note that a similar conclusion is reached by Sales et al. 



the disk plane, when it finally settles into the disk it only 
keeps the vertical component of its angular momentum. 
It then settles into a smaller, more concentrated disk. 

It thus seems that the formation of bulgeless galaxies 
requires an extremely calm gas accretion history, both 
in term of accretion rate and in terms of changes in 
the direction of the accreted angular momentum. The 
strongest correlation seems to be between the bulge con- 
tent at z = and the gas accretion rate at high redshift. 

To investigate the reasons why a more intense and more 
bursty accretion history gives birth to bulge-dominated 
galaxies, we study in more detail the five galaxies that 
have a high gas accretion rate at high redshift. Among 
these 5 galaxies, one has a major merger between z = 2 
and 1 so that gas accretion is probably not the main 
reason for the high bulge content. Among the four re- 
maining galaxies, a visual inspection of i-band images 
at z = 1 reveals that three are cases of very early bar 
formation, and the last one is undergoing violent disk 
instabilities, eventually resulting in the formation of a 
bulge and finally a small bar a few 100 Myr after the end 
of the clumpy phase (this galaxy is shown in Figure 22). 

Interestingly, face-on images of these 4 galaxies still 
show a bar at z = 0. However, three of them have a low 
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Fig. 22. — Disk instabilities in one of the galaxies with a high gas accretion rate at z > 1. The images show B-band surface brightness 
maps (30x30 kpc, color scale from 18 to 25 mag arcsec - 2 ) that are 500 Myrs apart, spanning a redshift range from 1.2 to 0.8. The initially 
smooth disk becomes unstable and fragmented, the fragments then migrate to the center of the galaxy and participate in bulge formation. 

value of e s (which means a low amount of rotation in the 
central regions), and these are the three galaxies with 
a Sersic index greater than 2. The fourth galaxy has a 
bulge with a lower Sersic index (n& = 1.1), but shows 
significant rotation in the central regions. 

The low values of e s could be due to a slow down of 
the bar with time (which h as already been shown in sim- 



o o. 



e o. 



ulations, see for in stance |Combes fc Elmegreen| 1993 
Athanassoula 2003). This slow down (and maybe par- 
tial dissolution) seems to be linked with the formation 
of a bulge with a high Sersic index. This picture, al- 
though based on an extremely small nu mber of cases, 
seem s in agreement with observations by |Durbala et aL 
(2008) showing that most of the local Sb-Sc galaxies with 
a bulge Sersic index greater than 1.7 have a bar. Our sim- 
ulations further suggest a link between old bars, formed 
in gas rich disks at z ~ 1 and bulges with a high Sersic 
index. 

To summarize this Section on the contribution of merg- 
ers and gas accretion to the bulge growth of simulated 
spiral galaxies (16 galaxies with B/T< 0.3): 

• we find a correlation between the bulge fraction at 
z = and the merger history, both quantified by 
the fraction of baryons brought in by major and 
intermediate mergers and by the mass ratio of the 
largest merger undergone after z = 2, 

• we also find a correlation between the bulge frac- 
tion at z = and the gas accretion history, partic- 
ularly when we consider the gas accretion rate at 
z>\ 

• the most disk-dominated galaxies both have an ex- 
tremely quiet merger history (with in most cases 
only minor mergers after z = 2) and an extremely 
quiet gas accretion history (gas accreted at a low 
and constant rate, with an angular momentum vec- 
tor always in the same direction). By contrast, 
more violent merger or gas accretion histories give 
birth to galaxies with more prominent bulges (Fig- 



ure 



23 highlights this result) 



• the galaxies with the highest bulge Sersic index at 
z = are not those with many mergers but those 
with intense gas accretion at z = 1 and either early 
bar formation or other disk instabilities. 

7. DISCUSSION 

7.1. Spiral galaxy formation in simulations 

Spiral galaxies formed in modern cosmological simu- 
lations usually suffer from two main issues: a too high 
bulge fraction, and a too high stellar mass for a given 
halo mass. Both issues can be found to some extent in 
our simulations. 
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Fig. 23. — Fraction of baryons brought in by intermediate and 
major mergers as a function of the average gas accretion rate at 
z > 1, with the simulated galaxies color-coded as a function of their 
final B/T. We find that galaxies with a low B/T have in common 
both a quiet merger history (less than ~20% of their baryonic mass 
brought by mergers) and a low gas accretion rate at high redshift. 
Both a more violent merger history or a higher gas accretion rate 
give birth to more prominent bulges. 



7.1.1. The bulge content of simulated galaxies 

Zoom cosmological simulations have recently started 
to be able to produce spiral galaxies with realistic bulge 
fractions, reaching values of 0.2-0.3 (|Agertz et al.||2011 



|Guedes et al.||2011| |Brook et al.||20liap . However' since 



gathering large samples of such simulations is extremely 
hard, it is still unknown whether these simulations can 
reproduce the observed range of bulge fractions. 

In our case, even if we do form some galaxies with a 
very low bulge fraction, and even if most of them host 
pseudo-bulges instead of classical ones, the distribution 
of our bulge-to-total ratios is significantly offset from ob- 
served distributions: we do not form enough bulgeless 
galaxies. 

This could in some cases be a consequence of the re- 
simulation technique we use, and in particular of the 
way we populate dark matter haloes with galaxies. In- 
deed, we assume that each halo only hosts one galaxy, so 
that when a satellite is accreted, satellites of this satel- 
lite are neglected. If the incoming stellar mass was dis- 
tributed not only in the accreted satellite but also in 
smaller dwarfs around it, these dwarfs could be more 
easily stripped during the infall phase, and the overall 
impact on the main galactic disk would be weaker. Our 
technique thus maximizes the impact of mergers, and we 
thus probably slightly overestimate bulge formation due 
to this process. However, many of the simulated galaxies 
have a quiet merger history, so that this effect alone is 
unlikely to account for the lack of bulgeless galaxies. 

This lack is possibly aggravated by the relatively low 
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number of violently unstable disks at z = 1-2 in the sim- 
ulations, which could instead lead us to underestimate 
bulge formation. Indeed, we find that ~25% of the sim- 
ulated galaxies show a ring morphology at some point 
between z = 1 and 2. These rings correspond to m=0 
instabilities, and in approximately half of the cases, the 
ring itself becomes unstable and fragments (quite simi- 
larly to what is shown in Figure 22), leading to bulge 



formation. In some other cases however, the rings do not 
fragment but dissolve slowly, and we also find cases of 
disks remaining very smooth and stable. If the stability 
of these disks were artificial, we would then underesti- 
mate bulge formation, and the lack of bulgeless galaxies 
would become even more stricking. 

These disks could be artificially too stable because of 
resolution effects. To test this, we have re-run one of the 
simulations at twice higher resolution (see the results in 
Appendix [A]) . In the standard resolution run, the chosen 
galaxy has a ring which remains stable for several Gyrs 
before slowly dissolving, leading to the formation of a 
bulgeless galaxy. At higher resolution, we find that the 
ring does not stay stable for such a long time, becomes 
unstable earlier, and a bulge grows more quickly. This 
gives birth however to only a very low mass bulge, with 
B/T ~ 0.02. We also find that the surface density profiles 
are very similar in the standard and high resolution runs. 

It thus seems that resolution alone could not be enough 
to explain the stability of some disks in our simulations. 
While it could also be due to our gas dynamics scheme, 
this seems similarly unlikely since such a scheme has al- 
ready been used successf ully to model gas rich unsta- 
ble disks (see for instance Bournaud et al.||2007a ). It is 
also probably not a consequence ol our implementation 
of supernova feed back (that could d issolve the clumps 



m some cases, see |Genel et al.||2012[ ), because our feed- 
back is rather weaE In turn, it might be linked with 
star formation, and with the overly rapid gas consump- 
tion that the galaxies seem to experience between z = 2 
and z = 1. At z = 1, the simulated galaxies only have 
gas fractions of 0.20-0.40, which might not be enough for 
strong instabilities to develop. 

A direct comparison with the observed fraction of 
clumpy high redshift disks is however hard to perform, 
particularly because of the different mass ranges ex- 
plored. A large fraction of our galaxies are small at 
z = 1-2, and have low stellar masses and SFRs (in the 
range 3-30 M /yr), which would make them hard to ob- 
serve at such high redshift s. 

7.1.2. The relation between stellar mass and halo mass 

In addition to forming galaxies with overly massive 
bulges, simulations also usually struggle to reproduce the 
observed relation between stellar mass and halo mass. 
Observational ly, direct measures can be obtained from 
weak lensing ( Mandelbaum e t al.| |2006 |van Uitert et al 



2011 ) or satellite kinematics (llvlore et al.|2uTT| . Many re 



cent studies have also used abundance matching, a more 
indirect technique where the observed stellar masses are 
associated with dark matter halos extracted from sim- 
ulations by assuming that the most massi ve galaxies 
are a s signed to the most massive (sub) halos (Guo et ah] 



2010} IMoster et al.||201Q| ). |Behroozi etaL] (2010) dis- 
cuss the various uncertainties associated with this tech- 
nique, arising mostly from errors in the stellar masses 



(mostly linked with stellar population modelling), with 
additional errors coming from the scatter in assigning 
galaxies to halos. Overall, most techniques g ive similar 
results within a factor of 2 (|More et al.||2011 ), with pos- 



sible differences as a function o f" galaxy Hubble type and 



color (iMandelbaum et al.||2006l Ivan Uitert et al. 



and with additional variations between central and satel- 



2011) 



lite galaxies (Neistein et al. 2011). It is interesting to 
note that at fixed stellar mass spiral galaxies appe ar to 
be living in lower mass halos compared to ellipticals ( |van| 
|Uitert et al.||2011[ ) . |Comer6n et al.| pOlT]) also point out 
that the stellar masses ol spirals could be underestimated 
by 10 to 50 % if the same stellar mass-to-light ratio is 
used for thin and thick disks, so that the stellar-to-halo 
mass ratio could be further decreased for these galaxies. 

However the stellar mass in simulations seems in most 
cases to be a few times too high (see fo r instance the 
compilation presented in Guo et aL]|2010 ). Our simula- 
tions suffer from the same shortcoming: depending on 
the simulated galaxy, the fraction of the cosmic baryons 
found in stars (i.e., the galaxy formation efficiency) varies 
between 0.38 and 0.88, with an average of 0.62. This is 
3 times higher than the observed fraction of ~ 0.2 (e.g., 
|Behroozi et al.|[2010[ ). Observational uncertainties prob- 
ably cannot account for such a discrepancy. Note how- 
ever that the agreement between simulations and obser- 
vations is usually improved when the stellar masses are 
measured in the simulations in the same way as observers 
do, i.e. when stellar masses are estimated from the s pec- 
tral energy distribution of th e simulated galaxies (|Oh| 
eTaT^[201l] |Guedes et aLpoTT] ) 



Overly high stellar masses could be due to the way we 
take the cosmic accretion of baryons into account, and 
could happen if gas accretion were overestimated. In- 
deed, our simulation technique assumes that filaments 
contain the cosmic fraction of baryons, which might not 
be true, and one of the consequences of the sticky par- 
ticle scheme is that we ignore hot halos and assume all 
gas is accreted cold. If less cold gas reached the galactic 
disk, this would help decrease the stellar content of the 
simulated galaxies. However, the galaxy formation effi- 
ciencies that we find are extremely similar to what sim- 
ulations using ot her gas mode lling technique s find (e.g. 
with AMR jAgertz et al.||2011l or w ith SPH, |Governato 
et al.| 2007| pcannapieco et al. 200 9} see also a co mpari- 
son of different codes in Scannapieco et al. 2012). This 
discrepancy suggest a much deeper issue, found m most 

simulations, whatever the technique used. 

An except ion has recently been published by |Guedes 



et al. (2011), whose disk galaxy does sit on the observed 
stellar-to-halo mass relation. However, they use a way of 
deriving stellar mass slightly different from what other 
simulations have used, and neglect metal- dependent cool- 
ing at high temperatures, thus probably maximizing the 
effect of supernovae. They argue that their success is due 
to a high star formation threshold, which they set at 5 
atoms cm -3 (this is 5 times higher than the threshold we 
use but it's unsure if this could be enough to cause the 
differences between the simulations). They present only 
one case however, and it would be interesting to see if 
their recipes also produce realistic galaxies over a wider 
range of masses and formation histories. Regardless, a 
high star formation thr eshold is probably not a u niversal 
solution by itself, since |Avila- Reese et al. (2011) used a 
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similarly high value but still find it hard to match the 
observed stellar to halo mass relation. 

7.1.3. Modelling the physics of baryons 

The fact that simulations struggle to reproduce the 
stellar masses and bulge fractions of observed galaxies 
suggests some deep issues with the current recipes used 
to model the physics of baryons. 

AGN feedback is probably not relevant for our study, 
given the halo mass range on which we focus, and given 
that black holes are probably not massive enough in late- 
type galaxies. In turn, star formation and supernova 
feedback parameters likely play a much more important 
role. While we already mentioned possible effects of the 
chosen star formation threshold, the star formation effi- 
ciency itself stro ngly impacts the m orphology of the sim- 
ulated galaxies ( |Agertz et al.|20lT ). Calibrating this effi- 
ciency is difficult, esp ecially since it cou ld vary with red- 
shift (as suggested by |Agertz et al.]2011 ), which would be 
the ca se if it depended on metallicity ( |Krumholz fc Dekel 
|2Q12| showed that taking into account the influence ol 
metallicity suppresses star formation, especially in small 
galaxies at high redshift). This calibration requires sim- 
ulations at much higher resolution than what can cur- 
rently be achieved in cosmological context, following the 
structure of the GMCs and capturing supersonic turbu- 
lence. Such simulations have indeed produced results 
differing significantly from lower resolution ones, showing 
for insta nce an increased star formation efficiency during 
mergers ([Teyssier et al | 2010|), or a different structure for 
the merger remnants ( Bournaud et al.|2011 ). 

Supernova feedback has been sho wn to be a key ingre- 
dient to produce realistic galaxies (iDekel & Silk||1986l 



and galaxy formation. 



Scannapi eco et al.||2QQ8l |Governato et al.| |2010||Fio 



ntek SteinmetzpQlip . Supernova- driven winds can- 
not escape the gravitational potential of massive galax- 
ies, which makes them inefficient at directly reducing the 
baryon fraction in Milky- Way type galaxies at low red- 
shift. However, some authors argue that the winds and 
galactic fountains are essential for redistributing angu- 
lar momentum and producing bulgeless galaxies, all the 
more so since most of t he ejected material has a low 



angular momentum (|Maller & Dekel||2002 
2Qllb|a 



Brook et al 

It seems unlikely that a stronger feedback 
would be the solution to the issues we encountered in 
our simulations, both because of the relatively low SFR, 
even at high redshift, in which case outflows are not very 
powerful, and because a significant fraction of the mass 
growth of the simulated galaxies happens at relatively 
low redshift (z < 1). 

In addition to these standard physical prescriptions, 
it could also be that current simulations are missing key 
physical ingredients like radiative transfer. Indeed, while 
radiation is usually thought to play an important role in 
regulating the fraction of baryons in low mass galaxies 
at high redshift (see for in stance the recent simulations 
by Petkova fc Springel|2010[ ), it might also be important 
at lower r edshif'ts and in massive galaxies. For instance, 



Cantalupo| (|2010j) argues that the cooling in a galactic 
halo can be regulated by the photoionizing radiation from 
the galaxy itself that would remove the main coolants 
from the hot gas, so that the cooling timescale would 
directly depend on the SFR of the host galaxy. This 
could be an efficient mechanism for regulating gas cooling 



7.2. 



The emergence of the Hubble sequence 



Observations have revealed significant differences be- 
tween galaxy populations at high redshift and the z = 
Hubble sequence. 

Already at z ~ 0.5 grand-design spirals are rare and 
the spiral structure of disk ga l axies is often more chaotic 

Disk galaxies are 
and at fixed stellar 



2002). 



van den Bergh et al. 2000 
still however observed up to z 
mass their siz e seems to undergo little change between 
z = 1 and ( Ravindranath et al.||2004{ |Barden et al.| 
2005). Large blue disks have also been found at z = 2 
(some with spiral arms), but these object are extremely 
rare, and probably corres pond only to the mo st massive 
galaxies at that redshift (iSzomoru 

Irregular and interacting galaxies become increasingly 
common with increasing redshift, with the population 
of such galaxies alrea dy significant a t z ~ 0.5. In the 



Hubble Deep Field, Abraham et al. (1996) find a frac- 
tion of irregular galaxies ol 40 % . ten times more than in 
the Local Universe. In the UDF Elmegreen et al. (2005) 
find only 31% of spirals and 13% ol ellipticals, the other 
galaxies being classified as chains, doubles, tadpoles or 
clump clusters. While most of these disturbed galaxies 
are probably disks (undergoing a phase of violent insta- 
bilities), they could not fit anywhere on the standard 
z = Hubble sequence. 

Most surveys conclude that the Hubble sequence dis- 
appears at z > 1.5-2, and most galaxies at that redshift 
seem to be irregular or compact ([Daddi et al. | 2004 Con- 
selice et al.||2005| |Papovich et al.pUU5p . 

Our simulations support this picture. At z = 2, the 
majority of our simulated galaxies have a stellar mass 
between 10 9 and 10 10 M . None of these galaxies has a 
spiral-like morphology: some have very thick disks, some 
are ellipticals and there is a large fraction of galaxies un- 
dergoing interactions. Their kinematics classifies them 
as dispersion-dominated systems and no clear Hubble se- 
quence could be defined. This may be a consequence of 
the mass range we explore, and is consistent with ob- 
servations showing that more massive galaxies are less 
irregular (with for ins tance some cases of m assive spirals 
as observed at = 2 by Szomoru et al.||2011 ). By z = 1, a 
number of our simulated galaxies have acquired a more 
standard spiral morphology, with a few cases of galaxies 
with spiral arms, and some with a bar. The majority 
of disks however show no signs of either a bar or spiral 
arms. Most galaxies contain several distinct dynamical 
components, with in some cases bulges and disks that 
can be separated from their kinematics. By z — 0.5 bars 
are much more common and the final structure is more or 
less in place even if spiral arms often appear fragmented. 
Between z = 0.5 and z — 0, spiral arms become bet- 
ter defined and bars continue their growth. Note that 
our simulations suggest bars to be long lived structures, 
which will be investigated in more detail in a future pa- 
per. 

We thus confirm that a sequence of disks and ellipticals 
is found at z — 1, but we find no one-to-one correspon- 
dence with z = 0. Most galaxies in our sample acquire 
their final morphology at z ~ 0.5 (or even later if the 
details of bars and spiral arms are taken into account). 
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7.3. Mergers and bulge formation in Milky Way type 

galaxies 

We notice that the most disk- dominated galaxies of 
our sample share extremely quiet merger histories, with 
only very minor mergers after z = 2. While this could 
seem quite obvious, it also contradicts recent claims that 
mergers could actually help disk survival by triggering 
central starbursts and subsequ ent strong outflows of low 
angular mo mentum material ( |Brook et al.|2011b Guedes 
etaLl[2QTT ) 



Our simulated galaxies show however that galaxies 
with a major merger at z = 1-2 can end up at z = 
with a relatively low bulge fraction of 0.15-0.3. In fact, 
we find that 6 out of 16 (i.e. 37%) of simulated galaxies 
with B/T< 0.3 have a major merger between z = 2 and 
z = 1. This f raction is significantly higher than the val- 
ues quoted in Weinzirl et al. (2009), where the different 
models tested conclude that less than 10% of galaxies 
with such a bulge fraction have a major merger at z < 4, 
which is a very strong constraint. Our simulations give 
much weaker constraints on the merger histories of spiral 
galaxies, provided that major mergers happen at z > 1. 

This result is a combination of some mergers having 
disky remnants (all of them have a Sersic index lower 
than 2) and of a significant amount of disk growth after 
z = 1. This significant disk growth at low redshift both 
comes from the accretion of fresh gas (although galaxies 
with mergers tend to accrete gas at relatively low rates — 
see middle panel in Figure [20t and from gas recy cling due 
to stellar mass loss ( |Martfgf"&; Bournaud||2010[ ). 

Regarding mergers and bulge formation, merger rem- 
nants can be significantly affected by the simulation res- 
olution : Bois et ah] ( |2010[ ) show that a resolution of 



~ 32 pc (compared to our resolution of 150 pc) is re 
quired to follow the evacuation of angular momentum 
necessary for the formation of slow rotators, and that 
the effect is even str onger for mergers of gas-rich galax- 
ies. |Bournaud et al.| ( |2011[ ) also showed that simulations 
resolving the turbulent and fragmented nature of the ISM 
produced merger remnants that were more compact and 
with higher Sersic indices. The low Sersic indices of the 
bulges formed in our simulations could thus be a conse- 
quence of their limited resolution. 

Independent to the detailed properties of the bulges, 
we find mergers to make a significant contribution to 
bulge formation at high redshift, all the more so since 
violent disk instabilities are relatively rare in our simu- 
lations. Whether this is realistic or not can be tested 
with studies of thick disks. Both the mass, shape and 
kinematics of thick disks are indeed tracers of their for- 
matio n mechanism ( |Bournaud et al. 2009| Sales et al. 
2009). The fact that t he edges ol thick dis ks are not 
fla red (see for instance Comeron et al. 2011) was used 
by |Bournaud et aL (2009) to show that they most likely 
did not form because ol minor mergers, but rather are 
the relics of gas-rich clumpy disks at high redshift. A 
future paper will be dedicated to thick disks in our sim- 
ulations, and if they are flared, this would indicate that 
we probably overestimate the importance of mergers for 
bulge formation at high redshift. 



8. CONCLUSION 

We present the results of a series of cosmological simu- 
lations targetting 33 galaxies in an isolated environment 
with halo masses between 2.7 x 10 11 and 2 x 10 12 M . 
We study the evolution from z = 2 to z = of their mor- 
phology with techniques based both on their kinematics 
and on 2D decomposition with GALFIT. 

We find at z = galaxies with a large range of Hubble 
types, from bulgeless disks to bulge-dominated galaxies, 
although the fraction of bulgeless galaxies may be lower 
than the observed fraction. Most of the galaxies host 
pseudo-bulges, with a Sersic index lower than 2, and 70% 
of them have a bar. 

At z = 2, the simulated galaxies are very perturbed, 
and if there is a disk, it is usually thick, and sometimes 
unstable and clumpy. By contrast, it is much easier to 
identify disks and spheroids at z = 1, even if spiral struc- 
ture is usually absent, and bars are rare (except for galax- 
ies with very high gas accretion rates at early times, in 
which case a bar can already have formed by z = 1). We 
find that spiral arms and bars are usually in place and 
are much more common by z = 0.5. 

Even if a Hubble Sequence could be defined at z = 1, 
we do not find any correlation between the morphology 
at z = 1 and at z = 0, with a whole range of possibilities 
for the z = 1 progenitors of spiral galaxies (interacting 
galaxies, bulge-dominated systems, pure disks, unstable 
disks...). There is a much better morphological correla- 
tion between z = 0.5 and 0. 

Focussing on the formation histories of galaxies with 
B/T< 0.3 (typically Sb and later types, and correspond- 
ing to about half the galaxies in our sample), we find 
a correlation between the bulge fraction at z — and 
the merger history, both in terms of the mass ratio of 
the largest merger undergone after z = 2 and in terms 
of the fraction of baryons brought in by major and in- 
termediate mergers. We also find a correlation with the 
gas accretion rate at z > 1. We note that the most disk- 
dominated of these galaxies both have an extremely quiet 
merger history (with in most cases only minor mergers 
after z = 2) and an extremely quiet gas accretion his- 
tory: they accrete their gas at a low and constant rate, 
with an angular momentum vector always in the same 
direction. 

By contrast, more violent merger or gas accretion his- 
tories give birth to galaxies with more prominent bulges. 
Interestingly, ~ 40% of the galaxies with B/T< 0.3 un- 
dergo a major merger between z = 2 and z = 1 (none 
of these galaxies has a major merger at z < 1). Their 
disk-dominated nature at z = is the consequence of 
both the relatively disky nature of the merger remnants 
and a significant disk growth at z < 1. 

Finally, we find that the galaxies with the highest bulge 
Sersic index at z = are not those with many mergers 
but those with intense gas accretion at z = 1 and either 
early bar formation or other disk instabilities. The rela- 
tion between old bars and high Sersic index bulges will 
be investigated in a future paper, but it se ems consistent 
with observations of local Sb-Sc galaxies by |Durbala et al. 
( 2008 ) showing that most galaxies with a bulge Sersic in- 
dex greater than 1.7 also have a bar. 

While our simulations are successful in producing 
galaxies with a large range of Hubble types, and even a 



19 



t = 4.7 Gyr 

O 

standard resolution 


t = 6.5 Gyr 


t = 8.4 Gyr 

• 


t = 10.2 Gyr 


o 

high resolution 


4 


* 


• 



10° 



standard resolution 

high resolution 



6 8 10 12 14 
Radius [kpc] 



Fig. 24. — Comparison of the stellar structure of a simulated galaxy in a standard run and at twice higher resolution. The left panels 
compare stellar surface density maps (30 x 30 kpc, identical colormap in all cases) for different times of the simulations, while the right 
panel shows the stellar surface density profile at t=10.2 Gyr. At both resolutions a ring structure first forms, but this structure is less 
stable at high resolution, and a bulge forms much earlier. In the end however, at both resolutions galaxies have a very similar structure, 
both in terms of visual morphology and stellar density profile. 

few nearly bulgeless galaxies, we still find that we prob- 
ably overestimate bulge formation as well as the stellar 
mass at fixed halo mass. These issues, that are common 
to many cosmological simulations, are very likely related 
to the modelling of the physics of baryons. A better cali- 
bration of the efficiency of star formation (with a possibl e 
dependency on metallicity, see Krumholz fc Dekel|2Q12 ), 
and a better treatment of stellar feedback (from sup er- 
novae explosions but also for instance the radiation pres- 
sure from young massive stars) are promising paths to 
follow in order to resolve these issues. In particular, the 
formation of disk-dominated galaxies seems to be easier 
in simulations where star formation is significantly de- 



layed ( |Scannapieco et aL 2Q12[ ). This effect should be 
convolved with our study to get a complete picture on 
how galaxies get their morphology. 
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APPENDIX 
RESOLUTION STUDY 

At high redshift, a number of our simulated galaxies show disks that remain extremely smooth, with some cases of 
rings that do not fr agment. Rings could be a frequent evolutionary stage in high redshift galaxies ( |Genzel et al.||2QQ8 



[Aumer et al.||2010| ), but they are usually thought to fragment into clumps, that can then migrate to the center anc 
form a bulge. Since we are interested in bulge formation, it is particularly important to make sure that our disks are 
not kept artificially too stable, in which case we would strongly underestimate the bulge content of our galaxies. 

The non-fragmentation of our rings is probably not an artifact of the sticky-par ticle technique, which h as already 
been successfully used to model gravitationally unstable disks (see for instance Bournaud et al. 2QQ7a| ). It could 
however be due to the limited resolution of our simulations. To test this, we have re-run one of the simulations with 
a twice higher spatial resolution (corresponding to 65 pc), and a mass resolution increased by a factor of 6 (the mass 
of a gas particle is then 2500 M ). 

We compare the structure of the simulated galaxy in the standard run and in the high resolution run. Figure [24 
shows that in both cases a stellar ring forms at high redshift. This ring is more unstable at higher resolution, it 
dissolves more quickly and a bulge is formed earlier. However, if the details of the time evolution are different, the 
final galaxies are very similar at both resolutions. At t=10.2 Gyr, both galaxies have acquired a spiral structure, they 
have similar radii and stellar masses: the total stellar mass within the inner 8 kpc is 2.17 x 10 10 M in the standard 
run, and 2.28 x 10 10 M at higher resolution. Furthermore, the surface density profiles are very similar, with the only 
difference being that the low resolution galaxy is bulgeless, while a small bulge has formed at high resolution. The 
bulge fraction at high resolution is however only 0.02. 

This test thus confirms that the stability of the rings can be an artifact of the resolution we standardly use. Increasing 
the resolution however only leads to a minor increase of the bulge fraction, so that most of our results should still be 
valid. 

THE SAMPLE OF SIMULATED GALAXIES 



In Figures 25 to 32 we present for each simulated galaxy: 



• i-band face-on images (70x70kpc, color scale from 16 to 27 mag arcsec -2 ) at z = 2, 1, 0.5 and 

• i-band surface brightness profiles at z = and the corresponding GALFIT bulge+bar+disk decomposition 

• stellar mass evolution with time (stellar mass within the optical radius) 
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Fig. 25. — Galaxies with a bulge fraction from 0.02 to 0.10 
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Fig. 26. — Galaxies with a bulge fraction from 0.10 to 0.19 
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Fig. 27. — Galaxies with a bulge fraction from 0.19 to 0.22 
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Fig. 28. — Galaxies with a bulge fraction from 0.22 to 0.30 
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Fig. 29. — Galaxies with a bulge fraction from 0.30 to 0.35 
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Fig. 30. — Galaxies with a bulge fraction from 0.35 to 0.50 
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Fig. 31. — Galaxies with a bulge fraction from 0.50 to 0.80 
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Fig. 32. — Galaxies for which no GALFIT decomposition has been achieved 
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